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Gaussian wavepackets are a popular tool for semiclassi- 
cal analyses of classically chaotic systems. We demonstrate 
that they are extremely powerful in the semiquantal analy- 
sis of such systems, too, where their dynamics can be recast 
in an extended potential formulation. We develop Gaussian 
semiquantal dynamics to provide a phase space formalism and 
construct a propagator with desirable qualities. We qualita- 
tively evaluate the behaviour of these semiquantal equations, 
and show that they reproduce the quantal behavior better 
than the standard Gaussian semiclassical dynamics. We also 
show that these semiclassical equations arise as truncations to 
semiquantal dynamics non-self-consistent in h. This enables 
us to introduce an extended semiclassical dynamics that re- 
tains the power of the Hamiltonian phase space formulation. 
Finally, we show how to obtain approximate eigenvalues and 
eigenfunctions in this formalism, and demonstrate with an 
example that this works well even for a classically strongly 
chaotic Hamiltonian. 



I. INTRODUCTION 

There has been a recent renewal of interest in the 
study of the quantum-classical correspondence problem 
in many fields of physics, termed 'postmodern quantum 
mechanics', in analogy with with the postmodern move- 
ments in the humanities |l| . These studies are an attempt 
to understand the characteristics of quantum systems in 
terms of the dynamics of associated classical systems 
If the classical system is obtained in the canoni- 
cal fashion by replacing operators with c numbers, the 
approach is the usual semiclassical approach. The as- 
sociation here is straightforward only when the classical 
system is integrable: the invariant torii of the classical 
phase space provide the Einstein-Brillouin-Keller (EBK) 
quantization rules 0. Matters are substantially more 
complicated when the classical system is non-integrable, 
which is the more general case. It is this Ti — > limit 
of the quantum mechanics of systems where the classical 
(h = 0) dynamics are chaotic that is not well-understood, 
and where the recent focus has largely centered . In 
particular, there has been work examining the statistical 
properties of the spectrum of such systems Q and the 
remarkable calculation method known as the Gutzwillcr 
Trace Formula, which says that the eigenvalues of the 



quantum system arise as poles in a weighted sum over 
the unstable periodic orbits of the classical system || . 

Chaos in classical mechanics is best observed in phase 
space, using techniques such as Poincare sections. Hence, 
alongside the spectrum, there have been studies which at- 
tempt to associate the complex structures of a classically 
chaotic phase space with the eigenfunctions of the quan- 
tum system. These have included the study of 'scars', 
which are the intensity peaks in an eigenfunction along 
the unstable periodic orbits of the classical system ||,[l0) . 
This idea also figures in the association of localized eigen- 
states with structures such as invariant stable and unsta- 
ble manifolds [|llj and cantori (BirkofT-Gustavson quan- 
tization (l|]). 

Neither quantum mechanics nor semiclassical mechan- 
ics has a standard phase space; all the analyses above 
must, perforce, carry out calculations in Hilbert space 
and then project onto the classical phase space using 
'symbols' such as the Wigner Distribution Function or 
the Husimi Distribution Function JTJ|. This makes a 
direct interpretation of the results, and the association 
of quantum properties with classical structures all that 
more difficult. 

There exists, however, another systematic approach 
to associating quantum mechanics with a classical sys- 
tem. This may be motivated in many ways, the simplest 
being the time-dependent variational principle (TDVP). 
This approach starts with the variational restriction of 
Schrodinger's equation to a subspace of the full Hilbert 
space. It can be shown that in this approach, one always 
obtains a 'classical' phase space, where the dimensions of 
the phase space correspond to the wavefunction param- 
eters. One can then derive and study the approximate 
(variationally guessed) dynamics of the wavefunction in 
a straightforward manner. Since this technique starts 
from the quantum Hamiltonian, with no reference to the 
classical system, and can be shown to exist even for sys- 
tems without a well-defined classical dynamics, it has 
been termed semiquantal dynamics fl4|||. In this paper 
we start with this semiquantal approach; we now sketch 
the motivation for developing it. 

Semiquantal dynamics studies a classical system (dif- 
ferent from the canonical one) to be able to understand 
the quantum system. Since these dynamics derive as 
a variationally restricted form of the Schrodinger equa- 
tion, it follows that as the size of the Hilbert space ac- 
cessed is increased, the full quantum dynamics will be 
more closely approximated by the classical system thus 
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obtained. [We point out that most numerical work in- 
volving 'exact' quantum mechanics uses a superficially 
similar argument, wherein the wavefunction is expanded 
in a basis set of finite size. There arc fundamental differ- 
ences between the two approaches; in particular, even 
when the Hilbert space is inherently non-compact, as 
in the systems we consider, basis set truncations are a 
finite-dimensional restriction to a compact subset of the 
Hilbert space whereas the semiquantal restriction is to 
a non-compact space.] There is hence intrinsic value in 
understanding these approximate dynamics. Further, as 
we will show, classical and semiclassical mechanics can 
be derived as special limits or truncations of semiquantal 
mechanics. This route thus provides a smooth transi- 
tion, in principle, between the two extremes of the clas- 
sical limit and the fully quantal system. The transition 
is made without, at any stage, sacrificing the access to 
the considerable formal power of a classical phase space. 
It therefore constitutes an alternative way to the under- 
standing of the quantum-classical correspondence ques- 
tion. This provides the formal motivation; since semiclas- 
sical dynamics is contained within semiquantal dynam- 
ics, the approach supplements the semiclassical proce- 
dure without being detached from it. Secondly, recognis- 
ing the connection between semiclassical and semiquantal 
dynamics enables us to propose a Hamiltonian semiclas- 
sical dynamics, which combines the advantages of both 
approaches. All the formalism of semiquantal dynamics 
carries over to this Hamiltonian semiclassical dynamics; 
this is our practical motivation 

Our discussion of semiquantal dynamics in this paper 
will be through the use of coherent states []T5|-|l7f . These 
arise naturally in semiclassical limits and also in the limit 
N — > oo where N is the number of degrees of freedom of 
a quantum system 18 Coherent states come in vari- 
ous flavors, and can be constructed in any quantum sys- 
tem if the dynamical group and the algebra are specified 
pFj|-p2"|]. Their simplest version is in ordinary quantum 
mechanics with a Hamiltonian H = p 2 /2 + V(x); they 
are Gaussian wavepackets. We work with this particu- 
larization; however, most of the ideas are constructively 
generalizable. 

Our presentation is organized as follows: in the next 
section we review the derivation of the semiquantal dy- 
namics of Gaussian wavepackets. In Section 3, we derive 
the standard semiclassical dynamics of Gaussians, fol- 
lowing the work of Heller |2^] . We then show explicitly 
that those dynamics, derived as consistent truncations 
in derivatives of the Hamiltonian, are not self- consis- 
tent truncations in the moments of the wavepackct. We 
thus propose an extended semiclassical dynamics 'mid- 
way' between these semiclassical and semiquantal dy- 
namics, which retains the Hamiltonian formulation, and 
hence constitutes a Hamiltonian semiclassical dynamics. 
We include a qualitative discussion of the comparative 
validity of these three dynamics in various regions of a 
classically chaotic phase space, and argue that the ap- 



proaches developed in this paper improve upon the dy- 
namics derived by Heller. In Section 4, we then develop 
the formalism of the phase space construction, including 
the explicit representation of arbitrary wavefunctions and 
of the propagator. Finally, we show how to obtain ap- 
proximate eigenvalues and eigenfunctions in this formal- 
ism, and demonstrate with an example that this works 
well even for a classically strongly chaotic Hamiltonian; 
we conclude with a short summary. 



II. SEMIQUANTAL DYNAMICS OF GAUSSIAN 
WAVEPACKETS 

The basic idea behind the use of Gaussian wave packets 
relies on the following I^J^pSf : 

1. Since Gaussians are their own Fourier transform, 
they are well behaved in phase space and their dy- 
namics has an invariant meaning; unlike the WKB 
method and the Van Vleck propagator || , there are 
no singularities that have to be dealt with by using 
'matching formulas' etc. 

2. Their dynamics are simply parametrized by the c 
number variables that specify the centroid (average 
variables) and spread (fluctuation variables). 

3. They form an overcompletc set; an arbitrary wave- 
function can be expanded in terms of Gaussians 
(albeit an infinite number of them) . 

This has lead to their use in methods known var- 
iously as the time-dependent Gaussian variational ap- 
proximation [p6|-p8|], the time-dependent Hartree-Fock- 
Boguliobov method or the squeezed coherent state 
ansatz. The usual derivation considers the time- 
dependent variational principle (TDVP) formulation 
p9|,p0[, wherein one posits the action 
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r = I dt< y,t\ih— - H\y,t> 



(2.1) 



The requirement that <5r = against independent varia- 
tions of < ^,t\, and \^f,t > yields the Schrodinger equa- 
tion and its complex conjugate (note that a ray rotation, 
i.e. \^,t >— > exp(i\(t) /h)\^ ,t > leaves the variational 
equations unchanged) . The true solution is approximated 
by restricting the choice of states to a subspace of the full 
Hilbert space and finding the path along which ST = 
within this subspace. 

While the formal derivation of semiquantal dynamics 
relies on the TDVP, we review here the alternative deriva- 
tion of semiquantal dynamics via the Ehrenfest theo- 
rem [ppljl; we believe that this is a simpler, somewhat 
more intuitive approach similar to formulations in non- 
equilibrium statistical mechanics. It also helps us make 
the connection with the standard semiclassical dynam- 
ics of Gaussian wavepackets in a straightforward fashion. 
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The dynamics obtained are the same as in the TDVP ap- 
proach. We work, for simplicity, with a one-dimensional 
problem. The extension to multiple dimensions is pre- 
sented in the appendix. 

Consider, therefore, a particle of unit mass moving in 
a bounded potential with a Hamiltonian + V(x) 

( V is assumed analytic in x and O denotes operators). 
Accordingly, the equations of motion for the centroid of 
a wavepacket are 



where 



— < x > = <p> 
at 



di <P> 



8V(x) 

< — > 

dx 



(2.2) 
(2.3) 



where the <> indicate expectation values. If V is at 
most quadratic in x then the <> factor; 



< 



dV(x) 
dx 



>= 



dV 

dx 



(2.4) 



and the centroid therefore follows the classical trajectory. 

The classical equations of motion do not hold for the 
centroid in the case of nonlinear gradients, however, and 
to study such systems we expand in Taylor series around 
< x > and < p > . The expansions are of the form 



< F(u) >= i- < Au n > F {n) 



n > 



(2.5) 



where F^ 71 ' — d n F/du n \ < u > , Au = u— < u > (we use 
the repeated index summation convention throughout 
this section of the paper unless otherwise specified). We 
note tha t if all the terms in the expansion are retained, 
Eq. 2J5 is an identity. Using this expansion and com- 
mutation rules of operators, we can generate a series of 
moment equations. In general these are infinite in num- 
ber; an arbitrary distribution (wavepacket) is completely 
specified only if all moments are known. We have little 
interest in the entire wave function and it would suffice 
to keep track only of < x > and < p >, thus creating an 
'Ehrenfest phase space'. The convergence properties of 
the infinite hierarchy are not well understood; it cannot 
be (or at least has not been) reduced in general . 

A possible way |3l| of obtaining a finite set of equations 
is to adopt a Hartree-like approach, requiring that all n- 
point functions be expressible in terms of the one- and 
two- point functions. This is equivalent to assuming that 
the wavepacket is a normalized squeezed coherent state 
of the form 



#(g,t) = (2?r/i) 3 exp 
and implies the relations 



l[A{q-xf +p-(q 



(2.6) 



(2mVu m 

< Ax 2m > = - — ^— (no summation), (2.7) 



m!2 m 
< Ax 2rn+1 > = 0, 
4/i < Ap 2 > = H 2 +a 2 , 
< AxAp + ApAx > = a, 



(2.8) 
(2.9) 
(2.10) 



A = — (i + a). 



(2.11) 



[This assumption is the same as that of the TDVP: the 
wavepacket is restricted to a given subspace. ] This yields 
the equations 



dx 

~db 

dp 

~dl 

dp 
~dt 
da 
~dt 



= P, 



i!2' 



7 V {2m+1) (x), m = 0,l,. 



a, 

h 2 



2/2 



(m- l)!2 m - 2 



V {2m) (a 



(2.12) 
(2.13) 
(2.14) 

to =1,2... (2.15) 



The system is now reduced to the dynamics of x,p,(i 
and a ( where we write x,p for < x >,< p >) and 
are exactly those derived from the action principle. To 
see the equivalence, equate the variables x,p,ji, and a 
wit h q, p, hG, and 4?iGTI respectively of the references 



|34|. that we do not, for the moment, have 

an equation for the phase of the wavepacket- this will be 
derived later.] 

If we now introd uce t he chan ge of variables /i = p 2 , 
and a = 2pII, Eqs. 2.12 through 2.15 transform to |3q] 



dx 
~dt 
dp 

~dt 

dp 

~dt 
dU 
~dt 



P, 



Jm 



-V (2m+1) (x), TO = 0,1,. 



m!2 J 

n, 

h 2 p 2 ™- 1 

V _ {m-l)\2 r - 



-V (2m) (x), to = 1,2. 



(2.16) 
(2.17) 
(2.18) 
(2.19) 



Remarkably, these new variables form an explicit gra- 
dient system, yielding an extended potential system as our 
approximation to the Hilbert space. The classical degrees 
of freedom are now the 'average' variables x,p and the 
'fluctuation' variables p, II, respectively; the associated 
Hamiltonian is 



Hp 



P_ 
2 



n 2 

2 



+ TT + Vext(x,p); 



(2.20) 



V ext (x,p) = V(x)+^ + ^—V^ 2m \x), to =1,2... (2.21) 



,2 m 



8p 



i!2 



where the subscript ext indicates the 'extended' potential 
and Hamiltonian. 

This formulation is very interesting and possibly quite 
powerful; the extended potential provides a simple visu- 
alization of the geometry of the semiquantal space. We 
may thus get a qualitative idea of the semiquantal dy- 
namics before proceeding to detailed analysis. 

We note here, for instance, that: 
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1. Both the fluctuation and average variables are 
treated on the same footing and the phase space 
is dimensionally consistent: p has the dimensions 
of length and II that of momentum, 

2. The value of H ext is < H > under this approxima- 
tion, and is conserved, 

3. V ex t has an infinite barrier at p — 0, arising, as ap- 
propriate, from the momentum squeezing. Hence, 
'quantum fluctuations' can not be zero except in 
the limit Ti — > 0: as stated, there are no caustics in 
this representation. 

4. V ex t is symmetric in p corresponding to a choice of 
sign in p = ^JJL = v 7 < X 2 >; the infinite barrier 
renders the choice one of convenience and of no 
physical significance. 

It is now straightforward to introduce a Poisson 
bracket on this space: 

r ^ i = dfdg | df dg df dg df dg 
' dx dp dp dH dp dx dH dp 

This bracket renders the dynamics into the standard form 



V = {y, Hext} 



where y represents the phase space variables x,p, p, U. 
The rate of change of an arbitrary function /(y) is 



dt 



— {f, H ext }. 



(2.24) 



This Hamiltonian picture for the evolution of the 
wavepacket parameters is not quite the complete one 
since, as before, a ray rotation leaves the expectation 
values unchanged. To determine the overall phase of the 
wavepacket, therefore, consider the Schrodinger equation 
itself 



where 



d 

dr 1 



|<I>,i >= exp(i\(t)/h)\V,t> 



(2.25) 



(2.26) 



If we substitute eq. |2.26| into eq. 2.25, and take the 
inner product with | 1 E r , t >, we get for X(t) the differential 
equation 



dA T . . d 
-=<*,t\ l H--H\%t>. 



(2.27) 



The right-hand-side of the above equation is the La- 
grangian corresponding to the extended Hamiltonian (see 
2.1). The phase is then just the action function T (di- 
vided by K) and its explicit time-derivative is: 



• pU-Ilp . 

A = \-px - H ext . 



(2.28) 



It is natural to separate the phase into two parts A 
A_d + Ag with the first part 



dr < H > 



-tH K 



(2.29) 



corresponding to the dynamical phase (i.e. a measure of 
the time of evolution). The second part is the geometrical 
phase: 



1 A S 9 

dr < in-— > 



dr ( 



dr 
pU - tip 



■px). 



(2.30) 



For the case of cyclic evolution of the wavepacket this 
corresponds to the 'Berry's phase' |36) of the evolution. 
It depends only on the geometry of the evolutionary path 
in phase space since it can be written in the canonical 
form 



X G (C) 



P dQ, 



(2.31) 



c 



(2.23) where P = ( p, n) and Q = (x, p). Note that the inte- 



grand of 2.3C differs from that of 2.31 by a total derivative 
which vanishes for the closed line integral. The expres- 
sion for Xg(C) has this simple form because of the Hamil- 
tonian character of semiquantal dynamics. We will use 
this later for the construction of eigenvalues and eigen- 
f unctions. 

The equation for the phase, along with the Hamil- 
tonian equations of motion for the evolution of the 
wavepacket parameters constitute the semiquantal dy- 
namics of the wavepacket. This lies on a space 1Z(2N) x 
5(1); for the case just considered, N = 2, the result is 
completely general, however. Finally, we note that the 
overall normalization factor for the Gaussian wavefunc- 
tion is a trivial function of p and does not consitute an 
independent variable. The normalization is maintained 
by the evolution. 

To summarize, the semiquantal restriction of the dy- 
namics to Gaussian wavepackets has taken us from a 
Hilbert space to a extended phase space and the dynam- 
ics are described on this phase space by a generalized 
Poisson bracket. An arbitrary operator of the Hilbert 
space may be associated with a function (called its 'sym- 
bol') through 



0(y) =<y\6\y> 



(2.32) 



where y = {x,p, p, li) and the \y > are normalized Gaus- 
sians serving as the basis set for the association. The 
dynamics of this symbol are determined by Eq. 2.24: 



the Heisenberg equations of motion are now represented 
by Poisson bracket equations. These are not the clas- 
sical equations of motion, however, since the generator 
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of time translations is not the classical Hamiltonian, but 
H ex t- We emphasize that for the expectation values of 
operators, the space TZ(2N) (the approximation to the 
projective Hilbert space) suffices for a complete descrip- 
tion of the dynamics. The dynamics of the wavefunction, 
however, also require the specification of the phase; this 
is what lives on K{2N) x S(l). 

As with any variational principle, these equations are 
approximations to the same extent as the restriction of 
the variational space. Similar equations hold for the full 
space: the actual specification of an appropriate varia- 
tional wavefunction may not, of course, be possible. 



III. COMPARISON WITH SEMICLASSICAL 
DYNAMICS 

In this section, we first derive the standard semiclas- 
sical dynamics of Gaussian wavepackets, following Heller 
|p3| (we hence term these dynamics, with apologies to 
others, Heller's semiclassical dynamics) . We demon- 
strate that though these equations arise as consistents 
truncations to scmiquantal dynamics in derivatives of the 
Hamiltonian, they are not self-consistent truncations in 
h. We thus are able to introduce an extended semiclassi- 
cal dynamics; this restores the consistency of the above 
truncation, and is a Hamiltonian semiclassical dynamics. 
We finally compare the three approaches - semiquantal 
dynamics, Heller's semiclassical dynamics and extended 
semiclassical dynamics. 

We make a prelimin ary n ote that any truncation to the 
sums in the equations 2.16 through 2 . 1S| at a specified or- 
der in m, say the nth, retains the Hamiltonian structure. 
Its properties (as developed in the next two sections) thus 
transfer over to the truncated systems. The only detail 
that one has to be careful about is the evaluation of the 
symbols for operators. This consideration and the inter- 
pretation of the truncations is simplified by the standard 
semiclassical ansatz that 



< x- 



9 9 

< x >> = p z 



o(h). 

This implies, therefore, that 

< p- <p» 2 = o{h) = n 2 . 



(3.1) 



(3.2) 



The truncations are then to the same order (the nth) 
in h. The calculation of the symbols is straightforward 
- one proceeds with the computation as if for the full 
semiquantal dynamics and then retains only terms to the 
specified order in h. 



A. Heller's semiclassical dynamics 

The semiclassical derivation of the equations of motion 
for Gaussian wavepackets, following Heller |2^| makes an 
expansion around the wavepacket centroid, truncated to 



second order in derivatives of the Hamiltonian. We make 
the same simplifications as before; the Hamiltonian is 

-2 

\ + V(x) and we use a one-dimensional wavepacket 



¥(g,t) =exp 



-[(<? - x)A(q -x)+p-(q-x)+. 



(3.3) 



The physical meaning of the parameters p(t), x(t) are, as 
before < x(t) >= x(t), < p(t) >= p(t) and A and s are, 
in general, complex. 

The fundamental approximation is a truncation of the 
Taylor series expansion of the Hamiltonian as 

H *-Y^f+ vi0) + v(1) ^ *) + \ y(2) ^ - ^ 

(3.4) 

where is the nth derivative of V(x) evaluated at 

x, as before. This is exact for quadratic potentials, of 
course. If we insert this, and the wavefunction above, 
into Schrodinger's equation, equate coefficients of (q — x) 
, and make the identification x = p, we recover classical 
equations of motion for the centroid: 



x=p, 

and obtain, for the wavepacket parameters 

02) 



A = -2A Z - 



ihA-V(x) 



El 

2 



px. 



(3.5) 
(3.6) 



(3.7) 
(3.8) 



At this point, using our knowledge of the derivation of 
semiquantal Gaussian dynamics, we make the transfor- 
mation 



s •■ 

A : 



s r 
a r 



ISi 

■ iai 



n 



n 



2h P V 

These then provide us the equations 



"V 



V(x) 



2 



px 





n 


Si = 


~p 


p = 


n 


n = 









(3.9) 



(3.10) 



(3.11) 

(3.12) 
(3.13) 
(3.14) 



along with 
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x = p, 

p 



-W 1 



(3.15) 
(3.16) 



as the equations of motion for the wavepacket. 

It is straightforward to show that Sj is the time- 
dependent normalization term [p3| . The solution for it is 
Si(t) = ^ln(p) modulo a constant; we ignore it hereafter, 
except to note that it is the same as the semiquantal 
case. The equation for s r can be rewritten, using the 
relationships 



V 2 >, 



pli = IT 



P n-Up 



where 



E 



2 



IP 

2 



V{x) 



n 2 

P_y(2) 



8p 2 



(3.17) 
(3.18) 

(3.19) 
(3.20) 



Note that this equation for the phase is the same as Eq. 
2.28 with E replacing H ext and that 



E = H„ 



(3.21) 



This system of equations has been used with great 
success to study semiclassical dynamics by Heller and 
co-workers [p3| and Littlejohn |25|, amongst others. In 
particular, it has been used to construct the semiclassical 
spectra of various sytems, by a method we shall sketch 
later. Littlejohn has argued Jl7j that the spectrum thus 
computed for integrable classical systems is exactly the 
EBK spectrum. For non-integrable systems, it provides 
an alternative to basis-set quantization and the Trace 
Formula. The interesting question of the connection be- 
tween these methods remains an open one. 



in to; there is hence no extended Hamiltonian for these 
equations and the formalism developed below does not 
apply. However, it is now easy to construct a semiclassi- 
cal Hamiltonian. 

We restore consistency in the truncation by adding 
back the 'extra' term ^-V^ 3 \ in the equation for p. This 
term arises naturally in the expansion around the cen- 
troid, but was neglected by the truncation to the second 
derivative. Hence, arguments similar to the above, trun- 
cated not in derivatives of the Hamiltonian, but moments 
of the Gaussian, give us an extended semiclassical Hamil- 
tonian 



H,. 



(3.22) 



This Hamiltonian is a truncation of the semiquantal 
Hamiltonian to m = 1 and by the arguments above, to 
0(h). This is what we call extended semiclassical dynam- 
ics - it lies 'midway' between the standard semiclassical 
dynamics and semiquantal dynamics. 

It is now simple to show that these systems contain 
within themselves the classical limit. If for each system, 
we take the physical limit p, n — > along with h — > 0, we 
get a Gaussian that corresponds in the limit to a classical 
point particle, and the dynamics reduce to that of the 
classical Hamiltonian. Equivalently, we could make the 
same 0(Ti) ansatz as before, and take h — > to recover 
the classical equations. 

We now compare the three approaches: semiquantal 
dynamics, Heller's dynamics and extended semiclassical 
dynamics and discuss their validity and utility. 



C. Similarity of the three aproaches 

We first wish to emphasize that all three approaches 
are similar in principle. The simplest way of seeing this 
is to consider their derivation in the Heisenberg picture, 
which we now sketch. This is a straightforward excercise 
that we do not belabor here. First rewrite the quantum 
Hamiltonian as follows: 



B. Extended Semiclassical dynamics 

How does this system compare with the semiquantal 
system we had deriv ed ea rlier ? It is a s imple excercise to 
show that equations 3.13 through 3.16 abo ve are trunca- 
tions to to = 1 of the semiquantal equations 2.16 through 



2.19 , except for a missing term ^-V^ in Eq. 3.16 for p. 
The phase, as mentioned, is also a truncation of the semi- 
quantal phase to m = 1. The structural correspondence 
betw een the two systems is then exact, except that Eq. 
3.16 is truncated to to = compared to the others . 



We have thus demonstrated that the standard equa- 
tions for the semiclassical dynamics of wavepackets cor- 
respond to truncations of the semiquantal equations de- 
rived earlier. These truncations arc to different values 



H= P - + V{x) 



■ pAp + Ap 2 + V(x) 



(3.23) 

V —Ax n V {n) (3.24) 



n=l 



=< u >,Au = u - u and = = 

<s> as usual. Now use this (still exact) form 



where u 
d n V/dx 7 

of the Hamiltonian to derive equations of motion for 
the Heisenberg operators x,p, Ax 2 and AxAp + ApAx 
(or any other operator for that matter). Having ob- 
tained these operator equations turn them into c-number 
equations by imposing a Gaussian ansatz and evaluating 
the expectation values of these equations in that ansatz. 
These equations are now the semiquantal equations de- 
rived earlier [the equations for the higher moments are 
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now redundant, as before, and may be safely ignored]. If, 
in the above expansion of the potential, we truncate at 
n = 2, and proceed as outlined, again imposing the Gaus- 
sian ansatz, we obtain Heller's equations. A truncation 
at n = 3 yields the extended semiclassical dynamics. All 
three approaches are hence identical in nature. In partic- 
ular, we note that all three are local expansions around 
the centroid. The global Hamiltonians H ext and H sc are 
effective Hamiltonians that the wavepacket sees, and are 
obtained by no more a global approximation than for the 
standard semiclassical case. The only point to consider is 
that the Gaussian ansatz seems 'natural' in the quadratic 
approximation; however, as far as the real system is con- 
cerned, it is just as valid (or invalid) a constraint for all 
the approaches. 

We also note that both the semiquantal and extended 
semiclassical dynamics can also be derived as examples 
of 'symbol' dynamics |37j]. In this approach, we project 
out a c-number symbol of the Hamiltonian (or its lo- 
cal quadratic approximation, respectively) in appropriate 
states (the squeezed coherent states), and treat this as a 
classical Hamiltonian, with an eye to understanding the 
quantum behavior. Heller's semiclassical dynamics are 
similar, but not Hamiltonian because the projection onto 
the squeezed coherent states is done after the dynamics 
are derived, instead of before, leading to an inconsistency 
in the truncation. 

For the semiclassical approaches, a time-scale of valid- 
ity [p8[ for the ansatz may be estimated and is normally 
argued to be C(Zn(?i _1 )), the 'log time'. Recent work 
jij, however, argues that this may be unduly pessimistic. 
Further, as we demonstrate below, many qualitative fea- 
tures of quantum behavior can certatinly be reproduced 
by these dynamics. 

D. Qualitative behavior in model potentials 

We now discuss how the qualitative behavior of these 
approaches conforms to intuition. We examine these dy- 
namics in a few simple model Hamiltonian correspond- 
ing to local neighborhoods of a classically chaotic phase 
space. We start with the corresponding quantal Hamil- 
tonian and derive from it the various dynamics. We then 
consider the behavior of wavepackets that start centered 
or slightly displaced from the origin in these model poten- 
tials. We will not go into the details; they can be easily 
worked out from the results of the last two sections. 

We first consider the more forgiving case of elliptic 
regions (equivalently, in general, an anharmonic well) 
where such a Gaussian ansatz might be expected to work. 
Our model Hamiltonian is 

H = ^-+x 2n (3.25) 

[where we consider n > 1 since all three dynamics are 
identical and exact in quadratic potentials]. The po- 
tential is bounded in x and its highest non-vanishing 



even derivative hence is positive. For the semiquan- 
tal dynamics, the form of the semiquantal Hamiltonian 
then guarantees that p always remains bounded as well; 
this is physically reassuring. However, this fails for the 
two semiclassical systems. In particular, for the sim- 
ple x potential, the dynamics for p for both semiclas- 
sical systems are unphysically unstable. This is easiest 
seen by considering the wavepacket centred at the origin 
(x(t) = x(0) — 0) where identically vanishes. The 
situation improves when the wavepacket is displaced from 
the origin though the dynamics are still asymptotically 
unstable with the extended semiclassical system perform- 
ing better than Heller's dynamics due to the coupling 
between the average and fluctuation variables. 

We now look at hyperbolic regions. If the model poten- 
tial is unbounded, the semiquantal approach will always 
give unstable solutions. This is not a breakdown of the 
method, but only as appropriate: an unbounded poten- 
tial must give a spreading wave-packet. However, note 
that for the unbounded Hamiltonian 

v 2 

H = ^+ax 2 -x 4 (3.26) 

(a > 0) both Heller's dynamics and the extended semi- 
classical dynamics, will be deceived into remaining stable 
with a centered wavepacket except in the singular limit 
a = 0. This persist even for small displacements with 
finite a for Heller's dynamics. For the extended semiclas- 
sical system the coupling of the average and fluctuation 
dynamics kicks in to drive the wavepacket away from the 
origin and give a spreading wave-packet. 

We turn to the more appropriate case of hyperbolicity 
embedded in a bounded potential, with a model Hamil- 
tonian hence, say, of the form 

H = ^--x 2 + x 4 . (3.27) 

If we once again start centered at the origin, the semi- 
quantal approach provides for a growing wavepacket that 
is ultimately reigned in by the quartic term, as is the 
appropriate physical behavior. Both semiclassical sys- 
tems will again be deceived - in this case into rapidly 
spreading without a bound. With small displacements, 
the semiquantal behavior is unaltered. Both semiclassi- 
cal dynamics improve marginally and can be stabilized, 
but in general remain unstable. 

These discussions demonstrate that in these model 
potentials, the semiquantal approach improves substan- 
tially upon Heller's semiclassical Gaussian ansatz, while 
the extended semiclassical approach does so at least 
marginally. The arguments above are not the whole 
story, of course, since the wavepacket will in general be 
moving through these regions of phase space instead of 
being statically centred there, but capture the essential 
features of the expected behavior. It should be clear 
that the coupling between the fluctuation and average 
variables is in general not simple to analyze. This is a 
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physically valid coupling and has its advantages, as just 
argued, but may lead to behavior which has to be inter- 
preted carefully. 



E. Interpretation 

In a recent paper pi, it was demonstrated that the 
semiquantal dynamics of the double-well system are 
chaotic, even though the classical and quantum problem 
are understood to be regular. This is to be viewed largely 
as a formal result: the transition between classical and 
quantum dynamics is not quite as straightforward as is 
usually understood; quantal effects do not always sup- 
press chaotic behavior. Further, an examination of the 
range of parameters and initial conditions within which 
chaos appears indicates the following realistic interpreta- 
tion. The dynamics of the wavepacket both deep within 
either well, and at energies such that the system looks es- 
sentially like an anharmonic well, are regular. It is in the 
transition regime, at energies close to the classical un- 
stable fixed point and separatrix, that the dynamics are 
complicated. It is the interplay between classically un- 
stable behaviour and quantal tunneling effects that are 
reflected in the observed chaotic tunneling. By the meth- 
ods outlined below, this would give rise to a rather com- 
plicated spectrum in this region - a wholly appropriate 
result. 

We reiterate that uncoupling the variables, thus yield- 
ing Heller's semiclassical dynamics, would do worse in re- 
flecting the quantal behavior in this system, as outlined 
above. Irrespective, it is obvious that this sort of behav- 
ior needs to be better understood. The discussion in the 
rest of the paper leads us towards doing so: matters do 
not end at having observed chaos or lack thereof in these 
classical-like dynamics. This structure is reflected in the 
spectrum and other pure quantal properties. 

Considering all these factors, therefore, we argue for 
the use of semiquantal and extended semiclassical dy- 
namics in exactly the same situations as the Heller's semi- 
classical Gaussian dynamics. Wavepacket dynamics must 
ultimately break down in the presence of chaotic classi- 
cal dynamics; however, the two Hamiltonian approaches 
should better reproduce the quantal behavior and thus 
improve upon the already remarkable results of Heller's 
semiclassical dynamics. The spirit of these approaches, 
we believe, is truly semiclassical; the Hamiltonian formu- 
lations make them particularly powerful. This is devel- 
oped in the following. 



IV. PHASE-SPACE PROPERTIES OF THE 
REPRESENTATION 

To understand formal power and properties of the 
Hamiltonian formulations constructed above, we briefly 



discuss the properties of coherent states using the stan- 
dard notation, simplifying our discussion considerably; 
we relate the notations later. 

If we consider the non-Hermitian operator 



1 



with 



a= -j={x + ip) 



a,a T = 1 



(4.1) 



(4.2) 



we can construct the states \a > defined as the eigen- 
states of a 



a\a >= a\a > 



(4.3) 



These states are minimum uncertainty Gaussians, known 
as coherent states. Due to the non-hermiticity of a, they 
are not orthonormal, nor are their eigenvalues real. The 
properties of coherent states [ ^o||3S|l are well known; we 
mention two we use later. Coherent states have the fol- 
lowing representation in terms of the number states, the 
energy eigenstates of the harmonic oscillator (in this case 
with m = h = u = I): 



\a >= 



,(- 



«f>. 



n > 



(4.4) 



and they form an overcomplete basis set. 

The coherent states do not include all Gaussians: 
we parametrize the space of all Gaussians (for a one- 
dimensional system) using four real parameters x,p, p, U 
whereas the coherent states use just two real parameters, 



(4.5) 



To understand the extra parameters, we consider 
squeezed coherent states. We follow Yuen in intro- 
ducing an operator b 



b = fia + va) 
with the c numbers [i, v satisfying 

H 2 -H 2 = i- 



(4.6) 



(4.7) 



The real and imaginary parts of this equation imply that 
fi and v correspond to two independent real parameters. 
It can be shown that 



[&,&+] = !, 



(4.8) 



and therefore that the transformation 4.6 is a linear 
canonical transformation. This provides b with prop- 
erties identical to those of a. In particular, there ex- 
ist a set of eigenstates of b ( which we denote \(3 > g = 
\a.\, a>2, fi, v > ) equivalent to the states \a >: these are 
the squeezed coherent states. These states include the co- 
herent states as a particular case, where fj, — 1 and v = 0, 
and in general include all Gaussian wavefunctions. The 
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independent parameters a±, Mi v i provide us with our 
four-parameter representation. 

If we consider the self-adjoint operators 0,1,02 defined 
through a = Si + W2 (corresponding to the quadrature 
components) we see that 



< Aaf >=< Ao, >= 



(4.9) 



(Am = to— < to > ) where all expectation values are 
taken with respect to the states \a >. Since [dj., d 2 ] = \, 
the uncertainty principle constraint on di,d2 is 



< Aa? >< Ao? > > — . 

16 



(4.10) 



The minimum uncertainty states \a > have equal uncer- 
tainty in the two quadrature components; their correla- 
tion is zero: 



< AaiAa 2 >= 0. 



(4.11) 



If we evaluate the same quantities in the states \ f3 > q , we 
get 



< Aa i ></ - 7lM-H 2 > 



< Aa| > q 



< AaiAa2 > q — —i{ix*v — v*n + 1) 



(4.12) 
(4.13) 
(4.14) 



The squeezed states, therefore, 'squeeze' one quadrature 
component at the expense of the other, and in general 
have correlations between the two components. How- 
ever, there is a linear canonical transformation of vari- 
ables (d — > 6) such that the new quadrature components 
b\ , 62 satisfy the minimu m unc ertainty equa tion in these 
states. We use here Eqs. 4.12 through 4.14 to make the 
connection with our previous notation, x and p are the 
same as above, and 



(4.15) 
(4.16) 



We now use these properties for the representation of 
arbitrary wavefunctions. 

We have so far considered the dynamics of a 
wavepacket that is constrained to remain Gaussian: this 
yields trajectories in a generalized phase space. However, 
it is also possible to construct propagators for arbitrary 
wavefunctions. In this case we have to consider functions 
on this phase space. We use the fact that it is po ssible 
to associate with every operator a symbol as in Eq. p2| 



and choose O to be the density matrix (or projection 
operator) P = >< \?|. We now associate with ev- 
ery wavefunction (upto an arbitrary global phase) a real 
positive function 



V{y) =< y\fy >< \fr|j/ > . 



(4.17) 



This function is an extended Q-representation (or ex- 
tended Husimi function) : we now demonstrate this, and 
touch upon its relevant properties. 

With the coherent states as defined above, the Q- 
function corresponding to a projection operator P is de- 
fined as |39|] 



Q(a,a*) = — < a\P\a > . 

TT 



(4.18) 



The Husimi function is defined in almost the same way, 
upto the normalization of the coherent states. Since the 
normalization of the coherent states is well defined, we 
use the terms interchangeably. 

The Q-function is essentially a probability distribution 
function, with useful properties for the computation of 
expectation values of operators, etc. [^9|. We do not 
use it directly, however, since the dynamics we have de- 
rived do not leave coherent states invariant. Hence, the 
dynamics in the Q-representation are awkward and the 
Hamiltonian formulation cannot be used to advantage. 
However, the space of squeezed coherent states is invari- 
ant under the derived dynamics. We take advantage of 
this by defining, as an intermediate tool, an extended Q- 
function over the space of squeezed coherent states. The 
semiquantal dynamics of this can be easily handled - we 
detail this below. The Q-function (and the projection 
operator itself) can also be extracted from the extended 
Q-function in a straightforward fashion - we first take 
care of this. 

We define the extended Q-function Q e xt{ct, p, v) as: 



Qext(a, A*) v) =< v\P\a, fi, v > 



(4.19) 



where we absorb normalization constants into the defini- 
tion of the states. We have also suppressed the depen- 
dence of Qext on the complex conjugates of a, /1, v. This 
function is precisely the representation we have derived 
through semiquantal dynamics: Q ex t = V. It is defined 
over a four parameter space a±, ct2, /x, v and has the obvi- 
ous and convenient property that a slice f44| ] of it, taken 
at \x = 4, v = 0, reduces to the Q-function: 



Q(a, a*). 



(4.20) 



This, and other expressions involving fi, v should be 
regarded as essentially formal ones indicating the various 
properties of the representation, rather than to be used 
for computational purposes. For that, we work directly 
in the x,p,p,Tl representation. 

We now come to the physical meaning of the extended 
Q-function: it is a generalized probability distribution 
function defined on this semiquantal phase space. To 
see this, we consider the properties of the Q-function, 
which we can obtai n as a slice of Q ex t- In particular, the 
definition Eq. 4.48 can be used to show that 
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d 2 a Q(a, a*) = 1 



(4.21) 



and that the averages of antinormally ordered products 
of the creation and destruction operators is 

<a r (a t ) s >= J d 2 aa r (a*) s Q(a,a*). (4.22) 

These are properties similar to those of classical prob- 
ability distribution functions. There are many discus- 
sions [jlj|,[39],[l4| on the merits of the Q representation 
versus those of others like the Wigner function and the P- 
representation. Each has been critically examined for its 
merits and demerits. We confine ourselves to the observa- 
tion that in most cases it can be easily shown that one can 
move between the various representations with the help 
of well-defined transformations, and that the representa- 
tions have properties which render each convenient for 
different calculations. We do, however, want to address 
the popular misconception that the Q-function (or the 
Husimi function) is an inadequate representation in that 
it cannot be used to reconstruct the original projection 
operator. As is known in the quantum optics literature, 
Q(a,a*) can be expressed as an absolutely convergent 
power series in a and a*. To wit, since Q(a, a*) can be 
written as 

Q(a, a*) = e (-°"-> ^ < nl ^-> a™(a*) n , (4.23) 
n,m f Vi!m! 

we can use the fact that for any projection operator P, 
the number state matrix elements satisfy 



< n\P\m > I < 1. 



(4.24) 



to show that the double power series for Q(a, a*) is ab- 
solutely convergent (see Gardiner |59| for details). 
If we write this power series as 



Q{a,a*) = ^Q„, m a m (a*) r \ 



(4.25) 



the number state matrix elements of P can be explicitly 
written as 



< n|P|m >= ir\/ n\m). Q n -r 



(4.26) 



We have also shown how to reconstruct the correspond- 
ing density matrix (or projection operator) from Q ex p- 
first project down onto Q{a, a*), and then use the power 
series expansion mentioned above. It should be pointed 
out that this is a consequence of the overcompleteness 
of Gaussians. A representation in a basis set with the 
usual properties of orthonormality would not permit this 
reconstruction. 

We also note that the eigenstates of x and p are con- 
tained as asymptotic limits of the squeezed coherent 
states, hence the reduced position and momentum dis- 
tribution functions of the wavefunction > may be ob- 
tained directly from Q ex t- To get this, we evaluate Q ex t 
along p — 5v with 8 real. The limit p — > oo and 5 —> 1 
is and p — ► oo, and 5^—1 yields |^(p)| 2 (see 

Yuen [B0[ for details) . This is simpler in the p, II repre- 
sentation. The two limits are now taken with II = and 
p — ► and p — > oo respectively. In essence, since Q ex t 
is defined on the space of all Gaussians, we have taken 
advantage of the fact that the ^-function can be obtained 
as an infinitely squeezed Gaussian [H. 

The extended Q-representation is hence a valid repre- 
sentation. Moreover, the the dynamics of arbitrary wave- 
functions in this representation are very simple. The von 
Neumann equation for the evolution of the density ma- 
trix is replaced by a Liouville equation, by exactly the 
same route that led to the Heiscnberg equations of mo- 
tion for arbitrary operators being replaced by Poisson 
bracket equations: 



Ot 



\$ext i Qext J 



(4.28) 



where H ex t is, as before, the expectation value of H taken 
with regard to the squee zed states. The Poisson bracket 



is defined by Eq. [2.22| . Thus, we have constructed a 
semiquantal and a semiclassical propagator: it is just the 
classical propagator for an arbitrary distribution function 



U(t) = exp{t{H ext ,}), 



(4.29) 



the exponentiation of the Liouvillian for the extended 
Hamitonian system. This propagator, by virtue of this 
construction, is well-behaved and has nice properties - in 
particular it satisfies the composition property 



u(h)u{t 2 ) = u{h + t 2 ). 



(4.30) 



Alternatively, P may be written as a normally ordered 
power series 



(4.27) 



Thus, the projection operator is obtained by replacing 
a m (a*) n in the power series for Q(a,a*) by the corre- 
sponding normally ordered product of the operators a 
and fit. 

This discussion therefore provides us with a way of con- 
structing an extended Q-function for any wavefunction. 



This is the primary formal result of this section. It is not 
quite a computational tool, however, since it corresponds 
to the calculation of an infinite number of trajectories. 

This discussion of phase space properties demonstrates 
the formal power of this Hamiltonian approach: In par- 
ticular, the construction of a propagator with desirable 
properties. In the following section, we show how the 
Hamiltonian properties help us to obtain eigenfunctions 
and eigenvalues. 
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V. QUANTIZATION 

We now sketch how to use these dynamics for requan- 
tization |43| in two different ways. The first, which uses 
the Fourier Transform (FT) of the survival probabilty 
(S(t) =< *(0)|*(t) >) is a computational method in 
that it is geared towards the numerical computation of 
eigenfunctions and eigenvalues, directly using the dynam- 
ics of the wavefunction. The second is a more formal 
consideration of the invariance properties of periodic or- 
bits and invariant torii, which enables the construction 
of a quantization rule associated with these structures. 
These separate considerations helps us bring out differ- 
ent properties of these dynamics. In particular, they help 
understand the spectrum obtained in terms of the phase 
space structures. 



A. Quantization using the survival probability 

As has been pointed out p3| , any approximate 
method of computing the time-dependent behaviour of 
a wavepacket can be used to construct a spectrum. The 
argument goes as follows. Note that an expansion of the 
wavepacket (at any time t) in terms of the eigenfunctions 
of the Hamiltonian exists, since they form a complete set: 



(5.1) 



All the time dependence is in the complex coefficients, 
since the states \<p n > are stationary states and the above 
can be written as 



Mt)>=Y, c ^°) e ~ iEnt/h \^> 



(5.2) 



E n is the eigenvalue of the appropriate eigenfunction : 

<<j> n \H\<j> n >=E n (5.3) 

This implies that the survival probability can be written 
as 

S(t) =< *(0)|*(t) >= K(0)\ 2 e- iEnt/R (5.4) 
by orthonormality : 



(5.5) 



It is now straightforward that the Fourier transform of 
S(t) yields the spectrum of the Hamiltonian. Extracting 
the eigenfunctions is equally straightforward; the Fourier 
integral of the wavefunction provides these: 



s«W/ft|*( t ) > =Cm |0 m > 



(5.6) 



The exact spectrum is available only when \^(t) > can 
be exactly computed. Any approximation can be used to 
generate a spectrum, however, which is in some sense 
as good (or as bad) as the approximate dynamics itself. 
When we use any of the extended dynamics above for 
the case of the harmonic oscillator, since these approxi- 
mations are exact, it is easy to show that this yields the 
exact spectrum and eigenfunctions [ [23|]3~ifl . For this case, 
the (x,p) variables decouple from the (p, II) variables, 
and all initial conditions give us periodic orbits (POs). 
We take advantage of these POs to evaluate (in this case, 
analytically) the infinite Fourier integrals above. [Note 
that this is already better than the WKB method, which 
does get the spectrum right (with a little bit of work to 
evaluate the Maslov indices), but does not yield the cor- 
rect eigenfunctions.] In general, for Heller's semiclassical 
dynamics, one has to work harder to get wave-packet 
quantization, even if one does find POs. This is because 
these POs are in the classical phase space and do not 
yield periodic behaviour in the fluctation variables. This 
is true even when they are stable POs. 

However for the Hamiltonian constructions, even for 
general systems, if a PO can be found in the extended 
phase space (irrespective of stability), the Fourier inte- 
grals for the spectrum and eigenfunctions can then be 
done exactly. This is a remarkable point: in principle, 
a single periodic orbit can be used to extract the en- 
tire spectrum of the system. There are some caveats, of 
course: a) Most POs are found through symmetry con- 
siderations; they cannot then provide the full spectrum 
and will only reflect the appropriate symmetry, b) The 
weight factors c m for most eigenfunctions will be small, 
for a general PO. Hence, it will not be practically possible 
to extract these eigenvalues and their eigenfunctions. In 
practice, therefore, several POs will be needed. There has 
been a lot of effort in recent years to find POs in Hamilto- 
nian systems, directed primarily at Trace Formula work 
p4j . This method is thus fortuitously situated. 

The features of the spectra obtained through the sur- 
vival probability method (working with Heller's semiclas- 
sical dynamics) are not fully explained however, [j45| . A 
primary problem is that there is has not been, so far, any 
means of verifying these calculations. The only compar- 
isons one could make is with spectra obtained through 
different methods (basis-set quantization, for example). 
The discrepancies could then be put down to the differ- 
ing nature of the approximations. Working within the 
Hamiltonian formulations of wavepacket dynamics, how- 
ever, enables us to construct another quantization pro- 
cedure: a different approach within the same approxi- 
mation. This second method is what we might term in- 
variance quantization. The mathematical arguments are 
essentially the same as those for Bohr-Sommerfeld quan- 
tization and are detailed elsewhere 
here a simple constructive version. 



46 1; we provide 



where we use the spectrum obtained from above to 
project out the rath state. 
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B. Berry's phase and quantization 

Consider the phase space trajectories obtained by solv- 
ing the ordinary differential equations (obtained above) 
for the parameters of the wavefunction. We can see read- 
ily that a periodic orbit solution to Hamilton's equations 
is invariant (on the 1Z(2N) subspace) under the action 
of the Hamiltonian. Thus, a wavefunction constructed 
as a sum over all the points of the PO (each point cor- 
responding to a Gaussian wave-packet) is a candidate to 
be an approximate eigenfunction for the system. This 
ignores the role of the phase factor: each point along 
the orbit acquires a phase factor during the evolution. 
The dynamical phase is no problem, since it is the same 
for all points along the orbit and can thus be factored 
as a global phase. The geometrical part Ag for an ar- 
bitrary PO, however, destroys the invariance in the full 
scmiquantal space 1Z(2N) x S(l). We must therefore con- 
sider the evolution of Ag as well. We now note that a 
PO x Ag such that the periodic evolution of Ag on S(l) 
is commensurate with that of the PO on 1Z(2N) consti- 
tutes a function invariant on 1Z(2N) x 5(1). This, then, is 
our eigenfunction. The requirement of the commensurate 
evolution of the phase translates simply to the relation 



P dQ = mh 



(5.8) 



\ G {PO) = I PdQ = n/i, 
Jpo 



(5.7) 



where we have used Eq. 2.31 and P = (p, II) and 
Q = (x, p) as before. This consideration of POs such 
that the point orbit acquires a Berry's phase equal to in- 
teger multiples of Planck's constant h has precisely the 
same form as the 'old' quantization rule of Bohr and Som- 
merfeld. However, since the rule applies in the extended 
phase space, as opposed to the classical phase space, it 
does not have the same meaning. In particular, there are 
no Maslov-Morse corrections M to this rule, since there 
are no singularities in the Gaussian representation. It 
can be shown pH , however, that the 'spread' variables 
p, n explicitly take care of these corrections. This had 
earlier been demonstrated through elegant group theoret- 
ical arguments by Littlejohn [[l7| for the case of Heller's 
dynamics. In general, the quantization condition above 
will give us results different from the EBK rule (the POs 
are in the extended space) but will always incorporate 
the Maslov correction. 

We can thus see how to quantize using these POs in 
the Hamiltonian approaches: the eigenfunction is con- 
structed as a weighted sum over the commensurate pe- 
riodic orbit. The weight factor for each point along the 
orbit is precisely the appropriate geometrical phase. The 
eigenvalue is H ext for that orbit (or H sc if we are using 
the truncated form). It is interesting to note that this 
quantization condition arises entirely from a considera- 
tion of Berry's phase for the dynamics. 

The extension of this argument from POs to invariant 
torii goes through easily 31 4(| and leads to a general 
quantization rule: 



where the closed integral is now taken over the ith irre- 
ducible contour around the torus and the quantum num- 
bers Ui are labeled accordingly. This is exactly Einstein's 
generalization [fi"8[ of the Bohr-Sommerfeld rule to invari- 
ant torii. It has been shown j4^^.^l[] that the eigenfunc- 
tions (or rather their Q-functions) are, by the construc- 
tion above, peaked around these invariant torii [ [SO] . 

We emphasize that this method can be used in either 
of the Hamiltonian formulations constructed above. We 
now outline an example ; the details are the subject of a 
forthcoming paper |31|. 



C. Quantization of a chaotic Hamiltonian 

As an example of the methods formulated above, we 
have studied the Hamiltonian 



H 



fy 



1 



(5.9) 



This Hamiltonian has recently been the subject of some 
large-scale numerical work |51|. These indicate that the 
classical version (O — > O for all operators) is a very 
strongly chaotic system, with few, if any, stable periodic 
orbits in the limit (3 — > 0. The quantal Hamiltonian with 
(3 = also resists numerical analysis: basis-set quanti- 
zation with matrices of dimension 3240 do not provide 
converging eigenvalues (see the paper by Eckhardt et al 
f3l| for details). They have hence used (3 = 0.01 for their 
analysis. As they point out, the eigenfunctions of this 
Hamiltonian belong to the symmetry classes of of the C± v 
symmetry group which has eight elements (4 reflections 
in the axes and diagonals and 4 rotations by 7r/2). The 
irreducible representations of this group split into four 
one-dimensional representation and one two-dimensional 
representation. They have restricted themselves to the 
four one-dimensional representations and used harmonic 
oscillator basis sets to obtain low-lying eigenvalues and 
eigenstates for this system. We have applied the semi- 
quantal method detailed above to this system. It works 
very well; in particular it does excellently for (3 = 0. The 
extended Hamiltonian in that case is 



H ext = -(pl+ P 2 y + n 2 x + nD 



(5.10) 



If we consider the same one-dimensional representa- 
tions as Eckhardt et al corresponding to wavefunctions 
which are A) symmetric under x — > y, x — > —x, B) an- 
tisymmetric, symmetric , C) symmetric, antisymmetric 
and D) antisymmetric, antisymmetric respectively, we 
find that all four subspaces can be studied by the sym- 
metry reduced version of H ext above: 
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(5.11) 



where (z,p) and (p, II) are the canonically conjugate 
pairs. This new Hamiltonian is explicitly integrable: it 
has two isolating integrals of motion, the total energy, 
and an invariant obtained by construction, 



expected quantal behavior and the extended semiclassi- 
cal dynamics do so marginally. Finally, we have shown 
how to obtain approximate eigenvalues and eigenfunc- 
tions in this formalism, and have demonstrated with the 
help of an example that it works well even for a classically 
chaotic Hamiltonian. 



(5.12) 



This fact enables us to find the invariant torii needed 
for the invariance quantization detailed above, and eval- 
uate the eigenvalues in a straightforward fashion. The 
details are provided in Rcf. |31|; however, we note here 
that the ground state energy we obtain analytically is 
Eq = 0.5953. This is substantially lower than the value 
obtained by Eckhardt et al of Efi um = 1.093. The semi- 
quantal method is based on a time-dependent variational 
principle; it is easy to show that for the ground state it 
reduces to the time-independent (Rayleigh-Ritz) varia- 
tional principle. The semiquantal result thus provides 
an upper bound to the exact result and hence improves 
upon the numerical results. This is true even when we 
use (3 = 0.01 in the semiquantal method. The ground 
state energy in that case is Eq = 0.6012, only marginally 
different from the (3 = case, and still better than the 
numerical result. 

This example shows that the semiquantal method can 
be used to advantage to obtain eigenvalues and eigenfunc- 
tions for classically chaotic systems and at least in some 
cases give excellent results. It also tells us what role of 
semiquantal dynamics is in the agenda of post-modern 
quantum mechanics: it is the search for a variational 
space such that the dynamics on this spa ce are explicitly 
integrable. We then apply the rule 5.8 to quantize the 
system, hence fulfilling the program of the extension of 
the 'old' quantization rule, and advancing the compre- 
hension of the quantum-classical correspondence limit. 



VI. SUMMARY 

In this paper, we have developed Gaussian semiquan- 
tal dynamics to provide a phase space formalism. We 
have thus been able to construct a propagator with de- 
sirable qualities. We have compared these equations to 
the standard semiclassical Gaussian dynamics and have 
shown that the latter arise as truncations to the semi- 
quantal equations not self- consistent in h, though con- 
sistent in derivatives of the Hamiltonian. This has en- 
abled us to extend the usual treatment and introduce 
a Hamiltonian semiclassical dynamics that retains the 
power of the phase space formulation. We have then qual- 
itatively evaluated the behaviour of these three different 
approaches in various situations and have shown that 
the semiquantal dynamics substantially improve upon 
the standard semiclassical equations in reproducing the 
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APPENDIX: GENERALIZATION TO N 
DIMENSIONS 

In this appendix, we discuss the generalization of the 
Hamiltonian formulation of semiquantal and semiclassi- 
cal dynamics to N dimensions. The simplest way of doing 
this is to use product spaces for the variational subspace, 
for example, 



(Al) 



With this factorization, all the one-dimensional argu- 
ments carry through exactly to N dimensions, and an ex- 
tended potential formulation (with a semiquantal phase 
space of 4 N) can be shown to exist |5l| . 

If we allow for correlations between the various de- 
grees of freedom, however, this is no longer possible. 
The Hamiltonian formulation still exists |]3C|| , but the ex- 
tended potential cannot be extracted. The best way to 
derive the equations of motion in this case is to use the 
TDVP and a wavepacket < x|W, t > of the form 



*(x,t) = exp 



■[(x - q) • A • (x- q) +p • (x- q) 



(A2) 



where now x, p, q are vectors and A a symmetric N x N 
matrix. If we continue to assume a classical Hamiltonian 
of the form 



Hdassical — TjP 2 + ^( x ) 



and rewrite A as 



A = i-G 



K 



(A3) 



(A4) 



with G, K real symmetric matrices, we can obtain the 
extended Hamiltonian 

H ext = H classical + Tr^G" 1 + K 2 G + GK 2 ) 

+ < V(x) > . (A5) 

where Tr indicates the trace operation on the matrices. 
The only calculation of any subtlety is the evaluation of 
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< V(x) >, since ^(x) is a multivariate Gaussian distri- 
bution. The result is a series, as for the one-dimensional 
case, in even derivatives of V. The coefficients of these 
terms are found by the standard prescription: 

< XiXjXk ■■■>=Y, < X " X b >< x cX d >■■ (A6) 

where the subscripts a, b,c,d,... are the same asi,j,k... 
taken two by two. The summation extends over all the 
ways in which i, j, k . . . can be subdivided in pairs. The 
terms on the right-hand-side of this equation are just the 
— £ N elements of G, which is the covariance matrix of 
the multivariate distribution. These elements constitute 
the canonical coordinates (along with the vector x ) of 
the extended Hamiltonian; the corresponding canonical 
momenta are the elements of K (and the vector p). If 
the two matrices are chosen diagonal, the system reduces 
to the product space discussed above, and the extended 
potential exists. It is the presence of the cross terms 
in K 2 G + GK 2 that prevents this in the general case. 
All arguments of the Hamiltonian formulation, includ- 
ing those about truncations and reduction to semiclas- 
sical dynamics carry through, irrespective: there are a 
few more equations, and they do not have as simple a 
form, however. We also note that the global phase of the 
wavefunction is, exactly as before, the extended action 
evaluated along the trajectory in the extended space. 
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